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' — " Stellar pulsations in rotating relativistic stars are reviewed. Slow rotation approximation 

is applied to solving the Einstein equations. The rotational effects on the non-axisymmetric 
■ oscillations are explicitly shown in the polar and axial modes. 
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. §1. Introduction 
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The amplitude of gravitational waves is quite tiny, although they carry enormous 
amount of energy. They propagate without suffering any disturbance after leaving 
the sources. The observation gives the direct information of the emission regions. 
The gravitational wave astronomy in the next century will therefore be an invaluable 
tool to diagnose the interior structure of the relativistic stars. As one of the promising 
| sources, stellar pulsations are important. The oscillations are excited in the newly 

born neutron stars just after supernova explosions, or in abrupt stellar quakes. The 
excitation mechanism is uncertain at present, much more the expected number of the 
sources. The direct observation only will give a clue to these questions. In addition 
to the dynamical oscillations, non-axisymmetric oscillations play an important role 
in rotating stars. They cause secular gravitational radiation-reaction instability, i.e., 
CFS-instability§) "H*. The instability was so far studied in the polar f-modes. See 
e.g., Ref. g). Recently, the same mechanism was examined in the axial r- modes W'H . 
The estimates based on the Newtonian calculations suggest that the instability sets 
in even for smaller angular velocity. The further refinement of the axial modes is 
^ ■ therefore crucial. The r-modes have been examined from various aspects™ ~c3 . The 

whole of the works can not be reviewed here. See the some papers and references 
therein for the details. 

In this paper, we consider the theoretical aspect of the stellar pulsations within 
general relativity Astrophysical applications of the CFS instability may be given 



elsewheretill. In §2, the mode classification of the stellar pulsations is given from the 
viewpoint of the physical forces and mathematical displacements. This is useful for 
the following discussion. In §3, the rotational effects on the polar modes are reviewed. 
In §4, a method is given to solve the axial modes in slowly rotating stars. In §5, 
the method is applied to the axial oscillations without metric perturbations. The 
approximation is known as the Cowling approximation E3* , and gives good estimates 
for some pulsation modes. In §6, the method is also applied to axial oscillations with 
metric perturbations. Finally, §7 devotes to the summary. We use the geometrical 
units of c = G = 1, in this paper. 



*' E-mail address: kojima@theo.phys.sci.hiroshima-u.ac.jp 
**' For a recent overview see the paper by Friedman in this volume 
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§2. Mode classification of pulsations 

2.1. Physical forces 

Pulsations reveal various involved physics in different modes of oscillation pat- 
terns and the frequencies. In order to explain the some pulsation modes, we consider 
linearized equations of motion in the Newtonian gravity: 

d t Svi + v k V k 5vi + 5v k V k Vi = -Vi (— ) - (ASP + t k A k V lP ) - V*<50, (2-1) 



P J P 

where r is adiabatic index, and A, is the Schwarzschild discriminant defined by 



P Fp 



(2-2) 



The first term in the right hand side of eq.(2-l) represents the pressure force. This 
driving force causes the p-mode. The frequency o~ p of the p-mode is typically given 
by the inverse of the free fall time, 

{ M\ 1/2 

*~{*) ■ (2 ' 3) 

where M and R are the mass and the radius of the star. The frequency increases 
with the number of nodes. The f-mode can be regarded as the lowest node-less one. 
The second term in the right hand side of eq.(|2-l|) is related with the buoyant force. 
The mode associated with the buoyancy is called g-mode. The typical frequency a g 
of the g-mode is characterized by the Brunt- Vaisala frequency uj bv > 



0~q ~ cjbv 



(-A-gf 2 , (2-4) 



where g is the local gravitational acceleration. The value A significantly depends on 
the position. The convectively unstable region corresponds to A ■ g > and stratified 
region corresponds to A ■ g < 0. For isentropic stellar model, we have A ■ g = and 
hence a g = 0. Since the thermal force is much weaker than the dynamical pressure 
in general, the motion driven by the buoyancy is therefore slow, and the frequency 
of the g-mode is much smaller than that of the p-mode, a g o~ p . 

These p-, g- modes exist in rather simple situation. If another physical ingredient 
is important, what happens? New effect changes the modes in general. For example, 
stellar rotation changes the frequency as 

a = a + o-'Q H , (2-5) 

where fl is angular velocity, and do is the frequency of the non-rotating star. In 
eq.( |2-5[) , we assumed the slow rotation and considered the first-order effect only. 
In addition to the shift of the frequency, the rotation causes a new driving force, 
i.e., Coriolis force. The new mode associated with the force is called r-mode. The 
frequency is given by 

a r ~ Q. (2-6) 
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In the limit of no rotation, the frequency is zero. In this way, a new physical effect 
causes a new mode associated with it. Another good example of the existence of the 
new mode is w-modeli-^ "Lr, which is related with the gravitational wave in general 
relativity. 

2.2. Oscillation pattern 

From the mathematical point of view, the oscillation modes can be divided 
into two classes by parity. The displacements on a sphere can be expressed by the 
'gradient' or 'rotation' of the spherical harmonics Y[ m . The polar (or poloidal) mode 
in the spherical star can be expressed by 

£ r = R(t, r)Y lm (0, </>), f = V(t, r)VY lm {9, 0). (2-7) 

The density and pressure perturbations are coupled in this mode. Another type of 
the displacements can be expressed by 

£ r = 0, i=U(t,r)(rxVY lm (d,4>)). (2-8) 

This class of displacements is called as the axial (or toroidal) mode. The density 
and pressure perturbations are zero on the spherical stars. Rotation however slightly 
induces them. 

These types of the displacements cause different types of the gravitational radi- 
ation. The change of the mass moment plays an important role in the polar modes, 
whereas that of the mass current plays a role in the axial modes. Both modes are 
subject to the CFS instability due to the radiation reaction. The polar f-modes work 
only in the rotating stars with nearly break-up speed. On the other hand, recently 
discovered r-modes work even in slow rotation. The axial mode therefore deserves 
to be examined more extensively. 

§3. Rotational effects on the polar modes 

In this section, the rotational effects on some p-modes are demonstrated. The 
formulation and the numerical results for normal frequencies of the slowly rotating 
stars are given in Refs. pH})- ^T| ). The non-radial oscillations with index I > 2 are 
described by the coupled wave-equations inside the star. They are schematically 
written as 

(-df + d*)X + F(d r X, d r Y, X, Y) = 0, (3-1) 

+ d 2 r )Y + G(8 r X, d r Y, X, Y) = 0, (3-2) 

where c s is the sound velocity. These two equations show that the perturbations 
propagate with the light velocity and sound velocity, respectively. The set of equa- 
tions is solved as an eigen-value problem with exp{— i(o~t — m(f>)}. We impose the 
out-going wave condition at infinity. The resultant eigen-value is a complex number. 
The imaginary part o~i represents the decay of the oscillations due to the gravitational 
radiation, if aj > 0. For the slowly rotating star, the frequency can be expanded with 
the rotational parameter e = Q^B? /M. As a result, the frequency is modified as 

o~ = <tr(1 + mfJ j? e ) — io~i(l + ma'je). (3-3) 
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Note that the axisymmetric mode (m = 0) is affected only from the second-order of 
the rotation. That is, the effect of the centrifugal force etc. 

Some frequencies of I = 2 mode are tabulated in Table |. The polytropic stellar 
model is adopted. The direct numerical results show that the corrections a' R and a'j 
are positive. This result means the counter-rotating mode m < beyond the critical 
velocity changes the sign in the pattern speed, i.e., the real part of the frequency and 
in the decay rate, i.e., the real part of it. This fact is the condition of the radiation 
reaction instability. In this way, the counter-rotating mode becomes unstable for 
large angular velocity e. Note that the f-mode is the most crucial among the p- 
modes, since the corrections are the largest and the f-mode becomes unstable first. 
The frequencies and the rotational corrections are calculated for a wide range of 
stellar models. The numerical calculations!!^ show that these corrections increase 
with the relativistic factor M/R. That is, fully relativistic calculation suggests that 
the critical angular velocity decreases. This means that the instability sets in even 
for smaller angular velocity, as the system becomes more relativistic. 

Table I. Normal frequencies of stellar model with n = 1,M/R = 0.2. 



mode 




/ 


a r R 4 /M 3 


o-/ 


f 


1.17 


0.57 


0.032 


3.10 


Pi 


2.70 


0.32 


0.006 


1.78 


P2 


4.12 


0.12 


0.001 


0.90 



§4. Perturbation scheme 

The axial oscillation is trivial in the non-rotating stars, since there is no restor- 
ing force in the fluid stars. The oscillation becomes possible in the rotating stars. 
We therefore have to consider the oscillation on the rotating stars, which will be 
a difficult task. We will consider the slow rotation approximation in the pulsation 
equations. That is, pulsation equations are expanded with respect to the rotation 
parameter. The method can be applied to the polar modes successfully. In the 
spherically symmetric case, the perturbations can be decoupled into the axial and 
polar perturbations with spherical harmonic index (l,m). They are respectively de- 
scribed by the axial functions Ai m = (Ui m , ho j m , hi ; m ), and the polar functions P[ m 
= {$Pim,Spim, Rim, Vim, Ho lm , H x i m ,H 2 i m ,K lm ). In the presence of rotation, the 
perturbations are described by the mixed state of them. If the perturbation equa- 
tions are expanded by the rotation parameter e, then the formal relation between 
the axial-led A[ m and the polar-led P\ m can be expressed as 

= [Am] + £ X [Piilm] + £ 2 X [Am, A±2m] +"', (44) 
= [Pl m ] +£X [Al±lm] + £ 2 X [P lmt Pl± 2m ] + • • • , (4-2) 

where the symbol £ denotes some functions of order e, and the square bracket for- 
mally represents the relation among perturbation functions therein. We assume that 
the axial-led and polar-led functions are expanded as 

A lm = A ( i i 2+e 2 A (2 2 + ---, Pim = + e 2 Pi 2 m ] + •••)• (4-3) 
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Substituting these functions into eqs.(4T)-(4-2), and comparing each order of e, 
we have the following equations of e n (n = 0, 1, 2), 

0=[42l, (4-4) 
0=[eV} 1 J lm + £xAi 1 2l (4-5) 
= [ £ 2 A®] + E x [eVg\ m ] + £ 2 x , A^J (4-6) 
= [ £ 2 A ( 2+^x^]. (4-7) 
We have here assumed that the perturbation is described by a single spherical har- 
monic in the lowest order, that is, Ap^ = 0, for /' ^ I, and used eq.(4-5) in eq.(4-6). 
Equation (4-4) represents the axial oscillation at the lowest order. Equation (4-7) is 

the second-order form of it, and the term £ 2 x A^ can be regarded as the rotational 
corrections. The method to solve the equations is straightforward. The first-order 
equations are solved by the axial-led functions. The polar-led functions are expressed 
using them. We have the second-order equations with the corrections expressed by 
the axial-led functions at the lowest-order. These equations are successively solved 
in the following sections. In the actual calculations, we also assume that the time 
variation of the oscillation is slow and proportional to ,!?, i.e., dt ~ fi ~ 0(e). This 
is true in the r-mode oscillation, as will be confirmed soon. 

§5. Axial oscillations in the Cowling approximation 

5.1. First- order solution 

We will apply the perturbation scheme described in the previous section to the 
problem of axial oscillations. In this section, we neglect the gravitational perturba- 
tion, that is, we consider the Cowling approximation. The details of the calculations 
are explained in Ref. [ll]). The leading order equation of 5T^ = is reduced to 

(&r-imx)ug>=0, (5-1) 

where 

x = WTtr = WTt) { °- uh (5 ' 2) 

and 8t denotes time derivative in the co-rotating frame, i.e., drUim = (dt+imf2)Ui m . 
It is easy to observe that the motion is trivial for the non-rotating case, i.e., dtU[ m = 
0. 

If one solves eq.( ^-l[) by the eigen- value problem, then the solution is expressed 
by a delta function. Instead, we will examine the evolution of the perturbation by 
the initial value problem. The solution of eq.( |5-l| ) is 

U^(t,r) = J fi 1 2{r) s+ .^ {Q _ x) ds = f^(r)e-^-^H(t), (5-3) 

where H (t) is the Heaviside step function. We will consider t > region only, so that 
the function H (t) may well be omitted from now on. The function describes the 
initial disturbance at t = . 
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The implications of the result fl5-3|) are as follows. In the Newtonian stars, the 
sinusoidal time-dependence can be described by a single frequency, 



^-x)-^=(l-^ I y) 



mQ. 



(5-4) 



This is the r-mode frequency measured in the non-rotating frame. In the relativistic 
stars, to is monotonically increasing function of r, wq < to < vjr. The possible 
frequency range is then spread out, 



1 



1(1 + 1) Q 



mQ < a < 1 



TO 



1(1 + 1) Q 



mQ. 



(5-5) 



In both cases, Newtonian and relativistic stars, the radial dependence fj' is arbitrary 

at this order. The function in eq.( ^-3[) is constrained by the equation of motions 
for the polar part, as will be shown in the subsequent subsections. In this meaning, 
the scheme (4-4)-(4-7) is degenerate perturbation scheme. 

5.2. Second- order equations 

According to (4-5), the polar-led functions (6pi±i m , 5pi±i m , Ri±i m , Vi±i m ) are 
expressed by U^. These corrections affect the axial parts with indices (7 + 2,m) 
and (l,m). We consider the equation with (l,m) only. Corresponding to eq.(4-7), 
the axial equation with the corrections up to 0(e s ) is 



= (d T -im X )U^+C[d T U ( ^}, 
where £ is the Sturm-Liouville differential operator defined by 

r 2 e (\-u)/2 



(5-6) 



8c3we 



(-A/2-1O 



(Po + Po) 



Av'(p +p ] 



(e^d T U^)' 



-(F + G)d T U^, 



(5-7) 



F 



-4c 3 w 2 e^ 2 -^ ^ 
4c 2 tu 2 e(- A - 3 ^/ 2 



rV/y y 
M ) 

e {\+3u)/2 



PO+PO 



Au'w 



(wr 2 e u )'(po+Po) 



+ 



(2c x e v 
\Av'r 2 



[(r 2 we 



-v\h2 



(5-8) 



G 



4C2 , 2 4 -u\l 



-3ci 
3m 



re 
2 



A/2 

6 



i(i + i) 



^6 

r 

Ito' 



3 to' 
2~to^ 



- + o— 6 + *2 
r 2 to 



5iy, 



e A 5W 3 

A; 2 -I m 2 H 

r to 



Wj | 6^3 \ 

TO TO / 



(5-9) 
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5.3. Solution of the radial function - mode specification - 

In order to solve eq, (|5-6|) , we introduce a complete set of functions y K (r) with 
eigenvalue — k with respect to the Sturm-Liouville operator £, 

C[y K }+Ky K = 0. (5-10) 
The second-order equation (|5-6| ) can be integrated with y K = ifnxfim-K( r ) an d un ~ 

(2) 

known function /y m (r) as 

Ug> = (im KX tfU K + f£>)e-+« a -&. (5-11) 
The sum of the first and second order forms is approximated as 

\l + imK X t)fW K + /«] e -im(n-x)t (5 . 12) 

'C + fl^e-^-^^, (5-13) 



u lm ^ u lm 



(2) 

where we have exploited the freedom of to eliminate the unphysical growing 
term in eq.(5-12). The value k originated from fixing of the initial data becomes 
evident for large t, since the accumulation of small effects from the higher order 
terms is no longer neglected. As a result, the frequency should be adjusted with 
the second-order correction to be a good approximation even for slightly large t, as 
eq.(5T3). This remedy is known as the renormalization of the frequency, or strained 
coordinate for t in the perturbation method©. 

We here summarize the calculations up to the second-order by the perturbation 
scheme (4-4)-(4-7). In the lowest-order calculation, the first-order spatial function is 
not determined, whereas the time-dependence is fixed. By considering the next order, 
the first-order spatial function fj^. K is specified with the second-order correction of 
the frequency k. The range and the nature of the correction n can not be explored 



without explicitly solving the eigen-value problem (5-10), which significantly depends 
on the equilibrium state, in .particular, A = or A ^ 0. For the barotropic case 
(A = 0), eq.flSIp is modified 0>. 



§6. Axial oscillations including gravitational perturbations 

6.1. Lowest-order calculation 

In this section, we incorporate the metric perturbations. The method to solve the 
axial oscillations is the same as in the Cowling approximation. We have additionally 
six metric functions (h lm , h x i m ,H i m , Hi i m ,H 2 i m , K [m ). In the lowest-order, the 
calculation is rather simple, since only two components, hp \ m and h\ i m are relevant. 
We define a function ^ m as 

^ = (64) 
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The relation between the metric functions is given by 



4 - 

re 



1 lm 



(8t — imw)$>\, 



lm 



(Z-l)(Z + 2 

The axial velocity function is expressed by two ways: 
(d T - imx) Uim = -47r(p + p)r 2 e~ u d T @i m , 



2imu' 

W+T) 



&lr 



u. 



lm 



Jr 4 v 



4^/ 
I in 



(v + l6ir(p + p)e x )<P, 



lm 



where 



J 



F [1(1 + 1) -2] 

-(A+i>)/2 



Eliminating Ui m in eqs.(6-3)-(6-4), we have the master equation as 
1 



(&r — imx) 



V&lr 



-Wirimxip + p)e ®i r 



This equation is reduced to an eigen- value equation l\ by 



assuming exp( 



{w-n) 



v<P 



lm 



Imi 



where 



/(/ + !) 
2m 



(a — mi?), 



q = — (jr 4 ^')' = 16?r(p + p)e x zu > 0. 



(6-2) 

(6-3) 
, (6-4) 

(6-5) 
(6-6) 

(6-7) 
-iat) : 
(6-8) 

(6-9) 
(6-10) 



Equation ( |6-8[ ) is called singular eigen-value equation, since it has a singular point tq 
unless q(ro) = 0, corresponding to the real value of p, = w(ro)- This kind of singular 
eigen-value equation is studied for the incompressible shear flow. See the Appendix 
for the Rayleigh equation in two dimensional parallel flows. Some important conclu- 
sion can be derived from the behavior of the background flow, i.e., the function w 
in this problem. For example, the necessary condition of the instability is that the 
function w has inflection point, i.e., q = 0. This condition is never satisfied inside 
the star. As a result, the frequency is real and the flow is stable in this order. 

6.2. From the Rayleigh to Orr-Sommerfeld equations 

In this section, the second-order corrections will be included. The calculation is 
straightforward, but very complicated in actual. We only consider the term with the 
highest rank of derivative with respect to r. The rank of the derivative is important 

*•* There is a misprint in the previous paper I- 1 . 



Gravitational Waves from Relativistic Stars 



9 



factor to determine the type of equations. The term is originated from the second- 
order derivative of Ui m as eq. (|5-6| ). Since XJ\ m is expressed by the second-order 
derivative of <Pi m , we have the fourth-order derivative of &i m . The term with the 
highest rank of derivative is explicitly given by 



imJ 



8c 3 - 



we 



(po + Po) r 
jAu 1 



2 r 



(Po+PoV' 



l in 



The second-order equation with this correction can be written as 



(dr — imx) 



(2)/ 

I'm 



-l6mm X (p+p)e x $W+V [<l>\ mi 



(6-11) 



(6-12) 



The term ( 6-11 ) effectively gives the 'viscosity' like the Orr-Sommerfeld equation in 
the incompressible shear flow. (See Appendix.) The viscosity is important for the 
stability of the flows. For the small Reynolds number, the laminar flow is realized, 
whereas the flow becomes turbulence above a critical Reynolds number. The effective 
Reynolds number R e in eq.( |6-12| ) is estimated from dimensional argument as 

R^^L^^L. (6-13) 

The viscosity term will play a key role on the singular point of the first-order equa- 
tion, but the consequence is not clear at moment. It is necessary to explore further 
how the effective Reynolds number should operate in the stability and so on. 



§7. Summary - pulsations and gravitational waves 



The polar modes, which exist in the spherical non-rotating stars, are extensively 
studied so far. The angular dependence of the modes is specified by a single spherical 
harmonic index. The decoupled radial equation can be calculated as an eigen-value 
problem. The rotational effects are also examined within the first-order, and the 
corrections are calculated as a = o-Q+me 1 ai + - ■ ■ . The relation between the pulsation 
and the gravitational radiation is evident in the modes, that is, the gravitational 
emission gives the imaginary part of the frequency. The gravitational waves would 
therefore give a good insight into the stellar interiors, when observed. 

Axial modes never exist in the non-rotating fluid stars, but exist in rotating 
stars. The oscillations are calculated for the slow-rotation approximation. When the 
mode is calculated by t 
are required in general 
dependence is dominated for a single spherical harmonic. As a result, the solution is 
constructed by an initial-value problem. The temporal dependence can be written 
as an infinite sum of the Fourier mode exp(— iai). The rotational corrections are 
calculated up to the third-order as a = e ai(r)+eras(r)+- • • . At present, the relation 
to the gravitational radiation is not clear. From the Newtonian estimates©, the 
radiation reaction affects the oscillations in order e 2m+2 . The reaction term should 
be even power of e, which is time-asymmetric and originated from the radiative 
boundary condition at infinity. Rather, the equations describing the axial oscillations 
are clearly related with the vortex as the Rayleigh or Orr-Sommerfeld equations. 



ie eigen-value problem, large number of spherical harmonics 
13. In this paper, we considered the mode whose angular 
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Appendix A 

Incompressible shear flow 



In this appendix, linear perturbation equations of the two-dimensional flows are 
summarized for the comparison. The detailed arguments can be written in some 
booksc^P of hydrodynamics. The basic parallel flow of incompressible matter is 
assumed to (0, u(x)). We consider the perturbation within the inviscid theory, and 
assume the stream function in the form as <P(x) exp{—ik(y — ct)}, then we have 



(u 



k 2 <P 



(Al) 



This is called Rayleigh's stability equation. Some theorems are proved for the equa- 
tion. For example, Rayleigh's inflexion-point theorem is stated as follows. A neces- 
sary condition for instability is that the basic velocity profile should have an inflexion 
point, i.e., u" = 0. However this condition is not sufficient. The sufficient condition 
is not yet found. 

Viscosity can not be neglected in some flows. If we consider the problem within 
the viscous theory, the perturbation equation of the two-dimensional flows is reduced 
to the Orr-Sommerfeld equation. Using appropriate normalization, we can express 
it as 

(d t + ud y ) (dl + = u "d y $ + R-\d 2 x + d 2 y f$. (A-2) 

In eq. (|A-2| ), R e is the Reynolds number defined by R e = uo£/v, where v is kinematical 
viscosity, and uq and I are typical scales of the velocity and length. Using the Fourier 
mode, we have 



{u 



k 



u"<P 



1 



ikR P 



dx 2 



k 



<f>. 



(A-3) 



Formally, the Rayleigh equation (Al) can be regarded as the leading term in the 
expansion of (kRe)^ 1 . The Rayleigh equation is a good approximation, if the second 
and fourth order derivatives can be neglected. When the condition is no longer valid 
in some regions, we have to take account of the viscosity term there. The stability 
of the viscous flow significantly depends on the Reynolds number. For the flow 
with low Reynolds number, the flow is laminar, whereas it becomes turbulent for 
the high Reynolds number. More detailed arguments with boundary conditions are 
indispensable to find out the critical number, flow pattern and so on. 
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